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Benign Paroxysmal Positional Vertigo (BPPV) is a mechanical disorder of the vestibular system in 
which calcite particles called otoconia interfere with the mechanical functioning of the fluid-filled 
semicircular canals normally used to sense rotation. Using hydrodynamic models, we examine the 
two mechanisms proposed by the medical community for BPPV: cupulolithiasis, in which otoconia 
attach directly to the cupula (a sensory membrane), and canalithiasis, in which otoconia settle 
through the canals and exert a fluid pressure across the cupula. We utilize known hydrodynamic 
calculations and make reasonable geometric and physical approximations to derive an expression 
for the transcupular pressure AP C exerted by a settling solid particle in canalithiasis. By tracking 
settling otoconia in a two-dimensional model geometry, the cupular volume displacement and 
associated eye response (nystagmus) can be calculated quantitatively. Several important features 
emerge: 1) A pressure amplification occurs as otoconia enter a narrowing duct; 2) An average- 
sized otoconium requires approximately five seconds to settle through the wide ampulla, where 
AP C is not amplified, which suggests a mechanism for the observed latency of BPPV; and 3) An 
average-sized otoconium beginning below the center of the cupula can cause a volumetric cupular 
displacement on the order of 30 pL, with nystagmus of order 2°/s, which is approximately the 
threshold for sensation. Larger cupular volume displacement and nystagmus could result from 
larger and/or multiple otoconia. 



I. INTRODUCTION 

Benign Paroxysmal Positional Vertigo (BPPV) is the most com monly diagnosed vertig o syndrome l)Brandtlll99l[K 
with a recent study suggesting that it affects 9% of older persons ( Qg halai et all |2000|) . BPPV is characterized by 
sudden attacks of dizziness and nausea triggered by changes in head orientation, and specifically afflicts the posterior 
canal. The disorder has earned the common name "top-shelf vertigo", since attacks often occur when the head is 
suddenly tilted back, such as when looking at objects on the top shelf. Other clinical features of BPPV include a 5-10 
second latency b etween the head til t and the onset of vertigo, and a fatiguable response which lessens with repeated 
head maneuvers l|Brandt fc Steddinl . il 993l) . Although the condition is not life-threatening, the disorientation brought 
on by attacks is severely discomforting, and can cause nausea, accidents and injuries. 

BPPV is caused by a mechanical dysfunction of the vestibular system in the inner ear, whose fluid-filled semicircular 
canals normally act to detect rotation via deflections of the sensory membraneous cupula. In BPPV, calcite particles 
(otoconia) are believed to interfere with the normal operation of the semicircular canals, falsely inducing a spinning 
sensation when in fact no rotary motion of the head is actually occurring. 

Two primary theories for the cause of BPPV have been advanced by the medical community: cupulolithiasis, in 
which otoconia are directly attached to the cupula, and canalithiasis, in which otoconia freely sediment through the 
canals and exert a fluid pressure on the cupula. A consensus is emerging that canalithiasis is the more likely mechanism 
for BPPV, supported in part by the largely successful clinical techniq ue used in BPPV treatment: therapeutic head 
maneuvers (B randt fc Daroffl 119801 lEplevi Il992 ISemont et all Il998[) designed to drive otoconia all the way around 
and out of the canal, so that the y settle in the vest i bule/ utricle. Further support for canalithiasis is reviewed in 
the work of Brandt and Steddin ijBrandt fc SteddirJ . [l9931. who compare the two mechanisms and conclude that 
canalithiasis is bett er able to explain the latency and fatiguability of the disorder. Recently, however, Buckingham 
(Buckingham, 1999) questioned the canalithiasis interpretation for the maneuvers, suggesting that settling particles 
should exert a transcupula r pressure from the vestibule as well as from the narrow duct. Finally, House and Honrubia 
(Hou se fc HonrubiaLl20f)lil) cite various clinical observations, some of which are consistent with canalithiasis, and some 
with cupulolithiasis, and suggest that in fact both mechanisms are viable and occur. 

Despite extensive quantitative modelling of normal vestibular functioning, the description and debate on mechanisms 
for BPPV have been purely qualitative. Without a quantitative analysis of these mechanisms, it is impossible to know 
whether either mechanism is p hysically capable of produ cing a response of the magnitude experienced in BPPV. 
Recently, House and Honrubia l)House fc HonrubiaL 12003^1 have taken an important first step toward a quantitative 
mathematical model of canalithiasis and cupulolithiasis. Specifically, they postulate that a settling sphere exerts a 
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FIG. 1 a) The human vestibular system consists of the semicircular canals for sensing rotation, the otolithic organs (saccule 
and utricle) for sensing linear/gravitational acceleration, and the cochlea for auditory sensing. Displaced otoconia, perhaps 
from the otolithic organs, are believed to collect in the posterior canal and cause BPPV. b ) Schematic of the model geometry 
for the semicircular canals (defined in Eq. Q), which is chosen to resemble measurements llCurthovs fc Omanlll987l) . A torus 
of major radius R consists of a thin circular duct of radius bd, and a thicker region consisting of the ampulla and the vestibule 
(located above the cupula, not shown here). The ampulla is spanned and separated by the membraneous cupula of radius b c , 
and the vector e|| points along the (local) canal centerline, from the ampullar side of the cupula towards the vestibular side. 
The position of a sedimenting otoconium is specified by its distance s p along the centerline of the canal and its distance r p 
from the centerline of the canal. The translating otoconium exerts a pressure field that displaces the cupula by a volume V c . 



transcupular pressure APh = ~F g ■ e||/37r6(s p ) 2 , where F g is the force on the particle due to gravity (corrected for 
buoyancy), e|| is the unit vector oriented along the canal axis, b(s) is the canal radius (which varies with centerline 
coordinate s), and s p is the particle location. This formula was argued based upon the following assumptions: (i) the 
component of the force perpendicular to the channel exerts no transcupular pressure; (ii) the fluid pressure is uniform 
across the canal cross-section containing the otoconium; (iii) the applied transcupular pressure arises only from the 
pressure-driven component of the sphere's drag, and not from viscous shear stress (whence the factor of 1/3); and 
(iv) (not explicitly stated) the Stokes settling velocity and pressure drag for a sphere are unaffected by canal walls. 

The present work builds on and extends the results of House and Honrubia in several ways. First, starting with 
the Stokes equations for viscous flow, we derive the transcupular pressure AP C resulting from a small particle set- 
tling through a fluid-filled channel. Significantly, both qualitative and quantitative corrections to APh exist, since 
assumptions (ii-iv) are not supported by solutions to the Stokes equations. 

Our expression for AP C follows directly from the equations of hydrodynamics, is independent of particle shape 
or even size (so long as the particle is small compared to the channel), and is valid for a canal of arbitrary (but 
slowly- varying) geometry. Furthermore, we consider two-dimensional motion of the otoconia (i.e. head maneuvers 
which are in the canal plane), so that otoconia settle away from the canal centerline. 

We simulate BPPV episodes by tracking the cupular volume displacement due to the transcupular pressure exerted 
by the settling otoconia. The cupular volume displacement is frequently assumed to be proportional to the 'sensed' 
rotation rate. Under such an assumption, the 'perceived' rotation rate (or, equivalently, the measured nystagmus) in 
our study can be obtained directly from the simulated cupular volume displacement. There is evidence, however, that 
additional neural processing ( called velocity storage) occurs in the b rain stem, and that the (measurable) nystagmus 
response reflects this process i|Raphan et all . Il979t iRobinsonl Il977|) . For completeness, we also simulate the nystag- 
mus response using a 'velocity storage' model of neural processing. We emphasize, however, that our fundamental 
picture and analysis of canalithiasis and cupulolithiasis do not depend on velocity storage. We compare our model of 
canalithiasis to a model of cupulolithiasis, and demonstrate that canalithiasis can actually give a stronger nystagmus 
response, with or without velocity storage. This agrees with the conclusions of House and Honrubia. 



3 



II. METHODS 

In this section, we derive a simplified model that includes the physics of settling otoconia and the processing of 
a velocity storage mechanism in order to capture the essential characteristics of canalithiasis, without requiring a 
detailed numerical investigation. Five components required for model simulations will be developed: (i) the basic 
functioning of the semicircular canals; (ii) an equation for the volumetric cupular displacement V c to an externally- 
applied transcupular pressure difference AP C ; (iii) the fluid pressure field set up by a settling particle of a given size, 
density, and position, giving the pressure difference AP C that distorts the cupula; (iv) the particle's settling velocity, 
accounting for Stokes drag, hydrodynamic interactions with channel walls, and background endolymph flow; and (v) 
the velocity storage mechanism. 



A. Overview of normal vestibular functioning 

The human vestibular system consists of two sets of organs: the otolithic organs, which sense linear acceleration 
and gravity, and the semicircular canals, which sense angular motion. Each ear has three semicircular canals that are 
oriented in a mutually orthogonal fashion in order to sense rotations about each of three axes (Fig. QJi). The canals 
are filled with endolymph of density p and viscosity \i similar to that of water. Each canal consists of a narrow duct 
connected to thicker ampulla and vestibule regions, which are in turn separated by the membraneous cupula. 

When the canal experiences an angular acceleration, inertia causes a lag in endolymph motion, distorting the 
cupula by a volume V c . The cupular distortion triggers activation of embedded sensory hair cells, signalling the 
angular acceleration/velocity to the brain and causing a compensatory eye motion called nystagmus. Under sustained 
rotation, the elastic restoring force of the cupula pushes endolymph back through the canal, and the cupula slowly 
relaxes to its undisplaced position. 

Neural signa l processing in the bra in stem modifies the transmitted signal in a process called velocity storage 
(Rap han et all Il979t iRobinsonl Il977j) . which lengthens the typical time constant of the associated sensation and 
nystagmus. The typical time constant for cupular relaxation (r c rj 4 s) is transformed into a nystagmus time constant 
of about t v w 21 s for rotation about the vertical axis of the head (yaw), where the lateral semicircular canal 
is stimulated. Velocity storage may be less effective in the posterior canal, or might also be less effective for an 
unnatural vestibular stimulation such as BPPV. Nevertheless, since experimental data regarding velocity storage in 
the posterior canal is presently unavailable, we assume it to be the same as in the lateral canal. 

The model geometry (Fig. QJd) we use in our simulations is chosen to closely resemble a semicircular canal — a 
two-part torus of major radius R, with a thick ampulla region and a thin duct region. We represent the inner and 
outer radii of the canal as functions of angle (0 < 9 < 2tt), 



Rout = R + b d and R in {0) = R 



A + B tanh 



(1) 



where A = 2b c + bd and B — 2b c — bd- The canal narrows from radius b c in the ampulla to bd in the duct roughly over 
an angle 25, centered around 9±. We define the cupul a to be located at 9 = 0, and choose numerical values (Table 
P) to reflect the measurements of Curthoys and Oman ijCurthovs fc OmanL Il987j) , with channel radii chosen to have 
cross-sectional areas equal to the roughly elliptical channels measured. We use (3dR and (3 U R (not shown) for the 
lengt hs of the duct and vestibule/utricle to be consistent with the geometry of Van Buskirk et al. ijVan Buskirk et al\ . 

The model geometry, though simplified, captures the most important shape variations from a hydrodynamic 
standpoint. 



B. Mathematical model for endolymph flow in the semicircular canals 

The classic mathematical mo del of the normal functioning of the semicircular canals, originally developed by 
Steinhausen l)Steinhausenl Il933|) . treats endolymph displacement with an overdamped pendulum equation. Two 
time scales appear in Steinhausen's equation: a rapid time scale 77, during which fluid moves ballistically through 
the canal, and a slow time scale t c over which the cupula relaxes back to its undistorted position. Subsequent 
theoretical work has typically worked towards calculating these time scales from the physical parameters and length 
scales characteristic of the vestibular system. Examples in clude the transient and q uasi-steady fluid motion in a 
torus with thin (duct) an d thick (vestibule /u tricle) sections ijVan Buskirk et allll976j) . ducts with non-constant and 
non-circular cross section l|Oman et a/.L [T987). allowances f or high-accelerations by cal culating the fluid flow near the 
cupula partitio n rather than as suming a Poiseuille profile l|Rabbitt fc Damiand . Il99~6|) . fluid flo w through a possibly 
porous cupula ijDamianol ^99) , and the effects of multiple hydrodynamically-connected canals ijMuller fc Verhagenl 
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Physical parameters for human semicircular canals 


R 


Major radius of semicircular canal 


3.2 mm 


fCurthovs & Oman. 1987) 


b d 


Duct radius 


0.16 mm 


11 


b c 


Mean radius of cupular partition 


0.68 mm 


11 


6i 


Angle where channel tapers 


36 ° 


11 


5 


Angle over which channel tapers 


5.4° 


11 


Pd 


Angle subtended by duct 


250° 


("Van Buskirk et al. 1976) 


f3u 


Angle subtended by vestibule/utricle 


75° 


11 


9 


Endolymph density 


1.0 g/cm 3 


fBronzino. 1995) 


P 


Endolymph viscosity 


0.01 g/cm/s 


11 


Po 


Otoconium density 


2.7 g/cm 3 


11 


a 


Otoconium radius 


0.5-15/im, 7.5 /j,m avg. 


fCamDos. 1990) 


Tc 


Cupular time constant 


4.2 s 


fDai et al. 1999) 


T v 


Nystagmus time constant 


21 s 


(Malcolm. 1968) 


K 


Cupular elastic constant 


13 GPa/m 3 


Calculated using r c in Eq. (|5J 


Ai 


First zero of Bessel Function Jo 


2.4 


("Van Buskirk et al. 1976) 


"th 


Threshold rotation for sensation 


2°/s 


("Oman et al. 1987) 


g 


Gravitational acceleration 


9.8 m/s 2 





TABLE I Physical parameters and derived quantities for BPPV. 



1988). Because the approach of Van Buskirk et al. gives results that are largely consistent with these revised models 
in the low frequency and velocity regimes considered here, we use a similarly simplified geometry in the present article. 
It should be straightforward to extend the present work to a more complicated (but realistic) three-dimensional model 
of the semicircular canals. 

Van Buskirk et al. explicitly calculated the time constants in terms of the system dimensions and parameters as 

where Ai is the first zero of the Bessel function Jq(X), i.e. Ai « 2.4. All physical parameters are known other than 
the elasti c constant K of t he cupula, here chosen to give the measured value (r c = 4.2 s) for the human cupular time 
constant l)Dai et all Il999(l . Following a step increase in angular velocity a(t) = f2o<5(i), the endolymph in the duct 
(and therefore the cupula) experiences a maximum displacement of volume 

_ 4pn o (l + /Vfl07ritta 

after a time Tf, after which the cupula relaxes back on a time scale r c . Using the values in Table the cupular 
volume displacement is related to rotation rate via £lo/V c w 5.6 x 10~ 2 °/pL s, and a threshold rotation rate of 2°/s 
corresponds to a cupular volume displacement of 35 pL. 

The influence of fluid inertia occurs over millisecond time scales, whereas the state of the fluid, cupula, and 
otoconium all change on significantly longer time scales. We therefore neglect fluid inertia throughout this study. 
With this approximation, the cupular membrane obeys a simple balance of pressures, in which three terms are 
important. A cupular membrane that is distorted to displace a volume V c exerts an elastic restoring force that 
contributes an amount —KV C to the transcupular pressure difference. Second, by conservation of mass, instantaneous 
cupular and endolymph volume displacements must be equal. The resulting viscous force contributes an amount 
—jV c to the transcupular pressure, where the viscous coefficient 7 will be calculated shortly. (Throughout this paper, 
dots denote time derivatives.) The last term to enter the pressure balance is an externally applied, time-dependent 
pressure difference AP c (i), which in our case is given by the fluid pressure set up by a settling otoconium. These 
three pressures balance, giving an equation for cupular volume displacement 

- KV C - Wc + AP c (t) = 0. (4) 

In normal vestibular functi oning, AP r arises fr om angular acceleration of the canals, in which case Eq. Q is consistent 
with Steinhausen's model l|Steinhausenl 11933^ if the effect of fluid inertia is neglected. 
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The viscous resistance coefficient 7 for quasi-steady fluid flow through the duct can be derived in a straightforward 
manner. Fluid in a straight channel of length PdR and circular cross-section of radius b^, subj ect to an a pplied 
pressure difference AP, moves with a parabolic velocity profile directed along the cylinder axis e|| IjLeal Il992~). 

uW = 3^« ('•<- <>''• < 5 » 

Note that we have assumed the usual no-slip condition for a solid/fluid interface. From the flow rate V c — 
irbjAP/8{i{3dR, the viscous coefficient 7 is shown to be 

7 = ^6^- (6) 

This result can be expected to provide a very good approximation for the low-Reynolds-number flow in the slightly 
curved canal duct and slowly-varying geometries of interest here. The pressure drop along the duct (which varies 
approximately with rj b^/ftd) is approximately a thousand times greater than that along the vestibule/utricle and 
ampulla (which varies with bf/[3 u ). Therefore, we approximate the viscous resistance in the canal as occuring in 
the duct alone. Furthermore, we approximate the channel as locally straight, because the radius of curvature of the 
channel R is large compared to the duct radius. Th is introduces errors of magnitude O (bd/R) ~ 0.05 in the axial flow 
and O (b d /R 2 ) ~ 10 -3 in the flow rate ijLeall Il99^ . Finally, although the channel radius is not constant, it typically 
varies slowly within the duct. 

The cupular volume displacement V c that results from an applied pressure AP c (t) is given by the solution of Eq. 

®, 

V c (t) = - ! AP c (Oe- (t -*' )/r ^', (7) 
7 J -00 

where r c is given in Eq. J5J. 



C. Pressure drop exerted by sedimenting otoconia 

The time-dependent pressure AP c (t) in Eq. JJJ) is provided by the stress field set up by a settling otoconium, which 
we model as a sphere of radius a <C bd and density p settling through the fluid due to gravity g. In Stokes flows, the 
velocity and pressure fields around a translating particle are insensitive to the detailed shape or size of the particle 
outside of the immediate vicinity of the particle. Rather, they depend only on the total force exerted by the particle. 
Therefore, we can approximate the fluid flow due to the particle as that of an equivalent 'point force,' since a < 6^/15 
even in the narrowest part of the curved duct. 

In Appendix^] we demonstrate that the difference in pressure between the fluid in front and behind a small particle 
settling in a (possibly curved) circular cylinder of local radius b(s) is given by 

AP C = 2%^i (l - J$) . (8) 
irb{s p y y b(s p ) 2 J 

Here en and s p are defined as in APh , and r p measures the radial distance of the particle from the canal center line 
(Fig. Qp). The gravitational force on the particle is given by 

F g = ^a 3 Apg, (9) 

where Ap = p a — p is the density di fference between the otoconium and the endolymph. Eq. (JHJ) was originally 
derived for a straight circular cylinder ijBrenneri 11958). but the physics of low- Reynolds- number flows allows it to be 
applied, as a first approximation, for channels whose geometry (e.g. radius or orientation) varies slowly. Due to the 
closed geometry of the semicircular canals, the pressure difference AP C set up by the flow due to the settling sphere 
is exerted across the cupula, so that the forcing AP c (t) in Eq. JJJ) is given by Eqs. JEJ and ©. 

The radial and axial positions (r p , s p ) of the otoconium determine the pressure AP C exerted on the cupula via Eq. 
(JHJ, to which the cupula responds via Eq. JJJ. Particle motion affects the pressure exerted on the cupula in three 
ways: (i) When otoconia fall from the ampulla (radius b(s p ) = b c ) into the narrow duct (radius b(s p ) — bd), the 
transcupular pressure is amplified by a factor b^/b 2 m 18; (ii) due to the duct's curvature, otoconia fall away from 
the centerline towards the walls (r p approaches b(s p )), which reduces the transcupular pressure due to the factor 
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1 — r 2 /b(s p ) 2 in Eq. Jjj}; and (iii) the unit vector e|| changes direction as an otoconium falls, which changes the force 
projection g • en. 

Eq. (JHJ is expected to hold with errors of order (a/b(s p )) 2 and b(s p )/R (up to roughly five percent), with three 
exceptions. First, the 'point force' approximation breaks down for particles located within a few particle radii of canal 
walls, which requires a more detailed analysis. The pressure is then expected to be small, but not zero as predicted 
by Eq. (JHJ). Second, the analysis leading to Eq. (JHJ) does not hold for particles near the cupula. In this case, the 
pressure distribution along the cross-sectional area of the cupula is not constant, and the total force on the cupula 
(i.e. pressure integrated over the cupular area) varies continuously from F g • e 1 1 (for a particle just below the cupula) 
to Eq. (JSJ) for a particle beyond a few channel radii from the cupula. We have not accounted for these variations in 
our simulations for small particles (Sec. IIII.Afl because most of the cupular volume displacemen t occurs whe n Eq. 
© is valid. Third, if the particle occupies a significant portion of the canal (e.g. "canalith jam" (Eplev, 1995)) the 
transcupular pressure would be greater than given in Eq. (JSJ). 



D. Sedimentation velocity of otoconia 

Here we give an approximate calculation of the sedimentation velocity of an otoconium, considering three possible 
corrections to the standard Stokes settling velocity: (i) particle inertia, (ii) advection with the background flow, and 
(iii) hydrodynamic interactions with the walls. The first two corrections are shown to be negligibly small, whereas 
hydrodynamic interactions are small except when otoconia are very close to channel walls. 

Since the particle is small co mpared to th e size of the channel, a good first approximation to its settling speed is 
given simply by Stokes' formula ijLealL Il99^1 

U. = = ^ « 0.2mm/s, (10) 

OTTfia 9/U 

using values from Tabic [I] 

The influence of particle inertia can be estimated by balancing the inertial term mx, where x(t) denotes position, 
against viscous resistance Qn^ax, from which an inertial time scale emerges 

to 2pa 2 in — 2 



67T/MI 9/i 



lCT^ms, (11) 



during which time a settling otoconium moves only a few nanometers. Thus inertia is negligibly small for settling 
otoconia, and is neglected here. Note that inertial effects such as the Magnus lift (on rotating particles) are likewise 
negligibly small. 

Next, we consider the influence of the background velocity of the endolymph, which is set into motion by the 
displaced cupula. A cupula displacing a volume per unit time V c sets up a parabolic flow in the duct (Eq. 10), 
whose maximum velocity is twice its average, or 2V C /Ttb 2 d . From simulations (e.g. Fig. |2J), we estimate a cupular 
volumetric displacement rate of order V c < 10 pL/s, giving a maximum endolymph velocity in the duct of order 
0.3/im/s — approximately a thousand times smaller than the Stokes settling velocity, given by Eq. (|10f> . We thus 
neglect endolymph flow in computing particle motion. 

Finally, we incorporate previously calculated hydrodynamic interactions between the particle and the walls of an 
enclosing cylinder. Details can be found in appendix [U] We list here the most significant features: 1) Hydrody- 
namic interactions are largely insignificant away from channel wall s, where the relative co rrection to Eq. (|10l) is of 
order O (a/b) < .08 unless the particle is very close to the wall (|Hirschfeld et all Il984j) . 2) Our approximations 
do not hold for otoconia v ery close to the wall, since non-spherical otoconia shapes and membraneous duct coatings 
l|Curthovs fc Omanlll987|) may play significant roles. Without reliable information about otoconium-wall interactions, 
we s imply assu me the otoconium to behave li ke a sphere, with near- wall parallel (U \\ ) and perpendicular (U\ ) veloci- 
ties ][BatcheloHll967tlCoolev fc O'Neill 1 19681 iGoTdman et all 119671 iHarroel fc Brenneii ll983UO'Neill fc StewartsonL 
119671) 

^11 ~ 7— StT and U±~U S -, (12) 
" log(a/a) a 

where d is the particle- wall gap. We introduce a cutoff d c = 1 /im for the particle- wall gap so that otoconia slide along 
the duct walls, as is required for the therapeutic maneuvers to be geometrically possible. Likewise, the transcupular 
pressure exerted by a particle very close to the wall is not known; however, we expect it to be small. 
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E. Velocity storage as a model of neural processing 

Our treatment thus far has concerned only the mechanical response of the system. Signal proce ssing is believed 
to occur in the brain stem in a process called velocity storage ijRaphan et al 1 119791 iRobinsonl Il977t) . While velocity 
storage does not affect the mechanical response of the semicircular canals, it is believed to affect the (measurable) 
nystagmus response, denoted E, and we therefore include it in the present work. (Another modification involves hair 
cell and afferent adaptation, which we do not treat due to its much longer (60-80 s) time scales). The effect of velocity 
storage is to lengthen the long time constan t of the horizontal semicircular can al t c « 4.2 seconds l|Dai et all 
to a longer time scale t v w 16 — 21 seconds ijDai et aZl Il999t fHaiii fc Zed. 11992ft . which improves the performance of 
the canals for angular rate sensing, extending their bandwidth to lower frequencies. The velocit y storage mechanism 
is typically desc ribed in Laplace transform space (indicated by tildes) using a transfer function ijRaphan et all Il979t 
IRobinsonl Il977| 

g(f) = 1 + stc (13 s 
f2(s) 1 + st v ' 

where fl is the angular velocity transduced by the canal (Eq. ©)■ Combining Eq. (|13fl with the Laplace transform 
of Eq. J3J, we arrive at an expression (in the time domain) for nystagmus as a function of transcupular pressure, 

^' ) ° 3mTft)g /l AJ> ° (t ' )e ''''''' /T, '' t ' (14) 

This expression allows us to translate the mechanical effect of a sedimenting particle into an equivalent sensation of 
motion, which is directly measurable as an eye movement. Again, however, we emphasize that the main affect of 
velocity storage is to lengthen the cupular relaxation time scale. The basic qualitative and quantitative conclusions 
we draw from our simulations are essentially unaffected by whether or not velocity storage is incorporated, other than 
the time constant for nystagmus decay. 



III. RESULTS 
A. Canalithiasis 

In this section, we describe simulations of canalithiasis in which we numerically integrate the motion of an otoconium 
settling under gravity. The orientation of the canal is held constant in these simulations, and the initial position of 
the otoconia is taken to be just below the cupula, since this represents the lowest point when the head is held upright. 
The transcupular pressure exerted as the otoconium falls (Eq. (jSJ) is calculated and used to determine the volumetric 
cupular displacement V c (Eq. Q) and nystagmus (Eq. ijlijll. We present here the results of a simulation of an 
average-sized otoconium starting just below the cupula in a model canal whose geometry (Fig. QJj) closely resembles 
the semicircular canals. In addition, in appendix iBl wc simulate a particle falling in a simpler geometry - a straight, 
tapering channel - to clearly illustrate the pressure amplification that occurs as the particle enters a narrowing canal. 

The trajectories of three otoconia (A, B, C) falling from the cupula into the curved tapering region (Fig. as 
well as the associated cupular volume displacement, pressure, and nystagmus (Fig. \Bp), highlight different features of 
BPPV attacks. Three distinct features can be seen in trajectory B: 1) During the first ~5 seconds, the otoconium falls 
through the ampulla and the cupular volume displacement grows only slowly; 2) During the following ^3 seconds, 
the otoconium falls into the narrow duct and the cupular volume displacement/nystagmus increases rapidly; and 3) 
Once the otoconium approaches the wall, it no longer exerts a significant pressure on the cupula, so that the cupular 
volume displacement decays with time constant r c . The peak nystagmus for a single average-sized otoconium is 
approximately the same with or without velocity storage: E max ss 2°/s when velocity storage is incorporated, which 
is to be compared with the purely mechanical response, AV c max « 30 pL, which also gives nearly 2°/s using Eq. 
Velocity storage mainly affects the time constant for nystagmus decay. 

Because we have imposed a minimum particle- wall separation d c , the particle moves significant distances even when 
close to the wall (Fig. We interpret the initial (« 5 seconds) period of low cupular volume displacement as 

'latency', and note that the delay in sensation and duration of latency predicted here are consistent with clinical 
observations. 

The trajectories of all simulated otoconia begin just below the cupula, because this location in the posterior canal 
is physically the lowest for normal head orientations, and is thus assumed to be an otoconium 'collection point'. No 
latency should be expected for otoconia started in the duct, although a significant nystagmus can result. 
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FIG. 2 Simulated canalithiasis episodes occuring when otoconia, starting at different radial positions just below the cupula, 
fall through a tapering toroidal channel, a) Otoconia trajectories, with dots to indicate the location of each otoconium at 
two-second intervals, b) Results for (i) the pressure exerted on the cupula for the three trajectories, (ii) the cupular volume 
displacement for each of the three trajectories, and (iii) the nystagmus response, reflecting the sensed rotation velocity. Note 
that Eq. can be used to convert AV C into nystagmus which gives peak values (£ max ~ 2°/s) that are approximately the 
same as those obtained with velocity storage. The largest cupular volume displacement occurs when the otoconium enters the 
narrow part of the canal, hydrodynamically amplifying the pressure it exerts by a factor of up to 36 times the naive estimate 
of Fg/nbc (which would occur during cupulolithiasis) . This occurs with a 6-13 second latency for trajectories A and B, and 
never occurs for trajectory C, which hits the wall before entering the duct, c) Simulations for particles of different radii, all 
started on trajectory B in the model semicircular canal geometry (Fig. Peak nystagmus and cupular volume displacement 
increase roughly linearly with particle size a for particles larger than a 3-5 microns. The time at which peak nystagmus occurs 
varies with 1 /a 2 . 



B. Influence of Otoconium Size upon cupular pressure and displacement 

Otoconia have been found in a range of sizes, typically between 0.5-15 /jm In this section, we 

provide some guidelines for the influence of otoconium size on the nystagmus and cupular volume displacement in 
canalithiasis. The results of simulations for different-sized particles (Fig. ^jp) demonstrate that peak nystagmus and 
cupular volume displacement increase linearly with particle size a for particles larger than a few microns, which can be 
understood as follows: Hydrodynamic interactions do little to an otoconium's trajectory until it comes very close to 
the wall, where the transcupular pressure is negligible. Thus particles of different sizes follow approximately the same 
trajectory, but with different speeds (U s oc a 2 ). Transcupular pressure is proportional to otoconium mass (AP C ~ a 3 ), 
but is only exerted while the otoconium is away from walls (a 'settling time' T s ~ l/U s oc a" 2 ). For particles large 
enough to settle faster than r c or t v , the exponentials in the integrals of Eq. Q or Eq. i|14|) . respectively, are 
approximately constant while AP C is appreciable. In this case, peak cupular volume displacement and nystagmus 
(oc T S AP C ) should increase linearly with particle size a (Fig. 



C. Cupulolithiasis 

In the cupulolithiasis mechanism, otoconia are physically attached to the cupula, causing it to distort when the 
head is turned, so that the pressure exerted by a head suddenly turned would be constant. The maximum cupular 
volume displacement and nystagmus would occur if the head were tilted and held in a position such that gravity acts 
normal to the cupula. In this 'worst-case' scenario, the transcupular pressure is given by AP C = (4ir/3)a 3 Apg/irb 2 ., 
and the corresponding cupular volume displacement is given by 

"•W-EBfv—-^)' (15) 

with a nystagmus response 



9 



IV. DISCUSSION 

Eq. ©, expressing the transcupular pressure set up by a small particle settling through a fluid-filled canal, is 
perhaps the central result of our work. It is significant that this expression follows directly from a solution of the Stokes 
equations for viscous flow. While House and Honrubia's APh does capture many of the important characteristics of 
canalithiasis, the form of that equation was argued based on a series of assumptions (some valid, some not) rather 
than from a solution to the equations of motion for the fluid. 

Eq. (JHJ) confirms several features of House and Honrubia's expression: (i) the characteristic magnitude of the 
transcupular pressure, (ii) the amplification of transcupular pressure when the particle enters a narrow section of the 
channel, and (iii) the inability of a force perpendicular to canal walls to result in a transcupular pressure. However, 
both quantitative and qualitative corrections to APh arise, as evident from Eq. JSJ: (i) the transcupular pressure 
AP C for a particle on the centerline is six times greater than given by APh - a surprising consequence of viscous 
resistance - and (ii) the transcupular pressure is diminished for particles located away from the canal centerline. 

Also significant is that Eq. © provides a quantative estimate (subject to the assumptions of small particles far from 
canal walls) for the transcupular pressure in an arbitrary (slowly-varying, gently-curved) geometry, without requiring 
a detailed numerical analysis of the resulting fluid flow. While we have employed a simplified model geometry, it would 
be straightforward to apply the present results to a more realistic, three-dimensional geometry of the semicircular 
canals, as well as to multi-step head maneuvers. 

House and Honrubia used APh in one-dimensional simulations with otoconia constrained to fall along the canal 
centerline, and found that an otoconium mass 0.087/xg (or 18 otoconia of radius 7.5/im) would be required to achieve 
a (measured) peak nystagmus of E w 42°/s. By comparison, an analogous one-dimensional simulation using Eq. 
(|5|) would require only such three otoconia (0.014 /xg) to achieve the same peak nystagmus. In two-dimensional 
simulations, however, simulated particles fall away from the centerline and diminish the transcupular pressure (Fig. 
|2J). Therefore, approximately 20 otoconia, or 0.1/zg, would be required to account for the peak nystagmus measured 
by House and Honrubia, since one average-sized otoconium can cause a nystagmus of approximately 2°/s. 

The observed latency of BPPV can be understood as the time during which otoconia settle through the ampulla 
(where the transcupular pressure is small), as was also suggested by House and Honrubia . As a result, particles that 
start in the duct rather than the ampulla should produce no latency l|House fc HonrubiaLl2003|) . 

During therapeutic maneuvers, the head is re-oriented several times so that otoconia travel around the duct and 
into the vestibule, a region common to all three canals (Fig. ^l). Two significant features can be expected: i) Towards 
the end of such a maneuver, otoconia traverse the crus commune, where the anterior and posterior canals are joined. 
During this period, the pressure field of the settling otoconia is exerted across both the anterior and posterior canals, 
and the direction of nystagmus should change correspondingly, from one that is in the plane of the posterior canal 
(upbeating mixed with torsion), to one in the plane of the sum of both canals (purely torsional), ii) As otoconia 
settle out of the narrow duct into wider regions (e.g. the vestibule or the crus commune), the transcupular pressure 
is reduced ac cording to Eq. (HQ). This pressure reduction provides a possible resolution to Buckingham's critique of 
canalithiasis ([Buckinghaml [l99^) : once in the vestibule, otoconia may still exert a fluid pressure on the cupula, the 
magnitude should be substantially decreased according to Eq. JSJ. 

Experimental investigations of BPPV models have used the v estibular systems of frogs. Otoconia were inserted into 
the open end of an isolated posterior semici rcular canal (PSC) llSuzuki et a/.l . ll996j) . or dislodged from the utricle and 
driven into the PSC of an intact labyrinth IfOtsuka et aZJ . 12003^) . Neural activity was measured following a change in 
PSC orientation, with free-floating otoconia ('canalithiasis') and attached otoconia ('cupulolithiasis'). For quanitative 
comparisons with the present theory, several ingredients would be required: physical parameters of the frog PSC, 
the relation between firing rate of PSC nerves and rotational motion, and the mass of otoconia in the experiment. 
Nevertheless, qualitative features agree with our model. In 'cupulolithiasis' experiments, little or no latency was 
evident, and the change in firing rate was sustained, with decremental time constants of 16.8±4.9 s and >43 s. In 
'canalithiasis' experiments, latency periods of 2-4 s were obserived, and decay time constants of 7-10 seconds were 
measured. 

I ■ . experimental results (|Suzuki et all Il996fl contrast with the predictions of the present th eory, the theory of 
House and Honrubia, and clinical observations ijHouse fe HonrubiaL 120031 ISemont et all ^998) in one significant 
aspect: a latent period was measured for a clump of otoconia started in the duct, which Suzuki et al. attributed to 
a re-orientation of the clump itself. It should also be noted that in these experiments, otoconia typically travelled in 
'clumps' and slid along side walls. 

The effects of particle inertia were estimated and demonstrated to be negligible from a hydrodynamic standpoint, 
which seems to run counter to the philosophy of clinical maneuvers designed with particle inertia in mind. Although 
inertia is irrelevant in the trajectories of small settling particles, it may play a significant role in dislodging otoconia 
adherent to canal walls. 

A large range of eye movement responses are possible for the various sizes, trajectories, and collective behavior of 
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otoconia. Multiple otoconia may be clumped or dispersed, and otoconia clumps may be treated by our methods so 
long as the clump is smaller than the duct radius. For example, a clumped group of N otoconia would behave like a 
single larger particle (of radius a' w A fl / 3 a), giving a maximum cupular volume displacement which is approximately 
iV 1 / 3 times larger than for a single otoconium. On the other hand, a group of N dispersed otoconia settling indepen- 
dently would result in approximately N times greater cupular volume displacement and nystagmus than for a single 
otoconium. 

Clumped otoconia thus cause less nystagmus than dispersed ones, a prediction that is seemi ngly at odds w ith 
the fatiguability of BPPV. Peak nystagmus typically decreases with repeated head maneuvers (Brandt, Il99l[) . a 
fact which is typically attri buted to the breaking and dispersion of particle clumps and/or particle margination 
l)Parnes fc Price- Jonesl [l993) . Several mechanisms might reconcile these pictures: i) Particle margination: repeated 
maneuvers presumably break up clumps and disperse the particles, which when isolated might adhere to the mem- 
braneous canal walls, ii) Exceptionally large particle clumps that occupy almost all of the duct would exert a larger 
transcupular pressure, and would move significantly more slowly, than individual particles. This would result in 
greater nystagmus, which would be lessened by the breakup of such a large clump, iii) The response of dispersed 
particles could result in weaker nystagmus than a single clumped stone if the dispersed particles follow different 
trajectories, since the nystagmus response is sensitive to particle trajectory (Fig. 

A comparison of the cupular volume displacement and nystagmus due to a single otoconium in both canalithiasis 
and cupulolithiasis suggests that canalithiasis is the stronger mechanism, and that multiple and/or larger particles 
are necessary to produce the same response in cupulolithiasis than in canalithiasis. While counter-intuitive, this 
follows from Eq. (JSJ for canalithiasis, which provides for a pressure amplification by a factor of up to 2b^/b 2 d ~ 36 over 
cupulolithiasis. In cupulolithiasis, if the head were tilted and held in place indefinitely, an average otoconium attached 
to the cupula would yield a cupular volume displacement V c ~ 2 pL and nystagmus response E w 0.6°/s, smaller 
than the corresponding canalithiasis results by factors of 10 and 3-4, respectively. Furthermore, the nystagmus due 
to a single, average otoconium in cupulithiasis is below the sensation threshold, meaning that multiple and/or larger 
otoconia are required for cupulolithiasis. 



APPENDIX A: An application of the Reciprocal Theorem for a point force in a tube of slowly varying radius 

In this appendix we use the reciprocal theorem from low-Reynolds-number hydrodynamics l)Leal ll992) to determine 
the pressure drop at l arge d istances from a point force (Stokeslet) in a tube. This analysis was apparently first given 
by Brenner ijBrennerl ^958) for a point force in a circular cylinder and we show that the same ideas can be extended 
naturally to a tube of arbitrarily varying cross section and arbitrary orientation relative to the direction of gravity. 

Recall that the reciprocal theorem relates any two solutions (u, a, f ) and (u, a, f ) to the Stokes equations 

Vu = (Al) 
V-cr + f = 0, (A2) 

where a is the stress tensor (with components Cy ), 

/ dui du-j \ ... 

^•=-^ + H^ + ^J- (A3) 

Here u and f represent, respectively, the velocity field and force per unit volume acting on the fluid. Then, within a 
volume V and corresponding bounding surface S, these fields are related by 

/ n • a • u dS - [ n • a ■ u dS = - [ u-fdV^+ / u • f dV , (A4) 
Js Js Jv Jv 

where n is the unit normal directed into the fluid domain from the surface. 

Here we consider the specific case where (u, a, f) represent the fields associated with a Stokeslet in an otherwise 
quiescent fluid inside an infinitely long tube of arbitrary cross-sectional shape, and where (u, a, f) represents the 
steady pressure-driven flow in a non-uniform duct. In this case, the volume V corresponds to the interior of the tube. 
Now, due to the no-slip condition on the walls of the duct, the surface integrals in equation (|A4|) are only nonzero at 
the 'ends' of the duct, which we consider to be at a very large distance (many cylinder radii) from the point force, 
where the Stokeslet velocity field u is negligibly small. 

We now simplify (|A4|> in steps. First, the pressure-driven flow contains no body forces in V, f = 0. Second, the 
flow due to a Stokeslet in a bounded cylindrical domain decays exponentially with axial distance from the singularity 
l|BlakeL 11979ft . so u — > towards the ends of the duct. Further, since f corresponds to a Stokeslet of strength F, then 
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f = FS(r), where <5(r) denotes the Dirac delta function, and we may evaluate the second integral on the right-hand 
side of HA4fl as u p (x p )-F, where u p (x p ) is the pressure-driven velocity field at the location x p of the Stokeslet. 

Since the 'ends' of the cylinder are considered far from the Stokeslet where u — > 0, then the corresponding stress 
field is simply a pressure that remains constant across each 'end' section. If we denote the difference in pressures 
between the two 'ends' as AP, then the first term on the left-hand side of (1) reduces to QAP, where Q is the volume 
flux due to the pressure-driven flow. Combining the above results gives 

AP = . (A5) 

Q 

This holds for a tube of arbitrarily varying cross-section and and arbitrarily-oriented force F. 

Finally, we invoke the lubrication approximation, which provides that for a circular tube of slowly varying radius 
b(s) the flow velocity is given by a generalized Poiseuille flow, 



fi = C/o(s) ( v 1 "6(^J e "' (A6) 

where s is the coordinate along the cylinder axis, and r is the radial coordinate, measured from the centerlinc. Using 
the corresponding flow rate, 

q _ nb(s) 2 U (s) ^ (A?) 
m Eq. {KB, we find the Stokeslet sets up a pressure difference 

AP.^A.^y. (A8 , 

Note that the exact location of the 'ends' does not matter, so long as the Stokeslet flow (u) is negligibly small at 
the ends. Eq. (JSJ follows from Eq. I|A8(1 by taking the ampullar and vestibular sides of the cupula as the 'ends'. This 
implies, however, that Eq. JHJ does not hold for particles 'close' (within a few canal radii) to the cupula, since in that 
case u differs from the assumed 'infinite tube' solution. For example, a Stokeslet loc ated much closer to the cupula 
than to the side walls exerts an average transcupular pressure F g /nb 2 ijBlakeL Il97l[) . which is half as large as that 
predicted by [HI 

Lastly, we examine the error introduced by the point-force approximation. Because variations in the canal radius 
are gradual, an error estimate can be obtained from results The pressure drop due to a sphere settling through an 
otherwise quiescent fluid in an infinitely-long, straight, circular cylinder, is given by Eq. l|A8jl . with errors of order 
O ((a /b) 2 )), when the sphere is not within a few radii of the wall (see Eqs. 7-3.96-7.3.98 in ( Ha ppel fc Brennerl 
1983)). This corresponds to less than one percent for an average otoconium. Corrections due to canal curvature of 
order bd/R, approximately 5%, may also be expected. 



APPENDIX B: Cupular response to a particle settling in a straight, tapering channel 

A straight, tapering channel provides a particulary simple illustration of the pressure amplification provided by 
a narrowing channel, which is expected to be a significant dynamical feature accompanying sedimentation of an 
otoconium, and highlights the hydrodynamic explanation for latency. 

We consider a particle of radius a that settles along the midline of a straight vertical channel of varying radius b(s) 

6(s) .^-fcMtanh(^), (Bl) 

so that b(— oo) = b c and b(oo) = bd (Fig. 1631 1). The specific functional form in Eq. I|B1|) has been chosen for 
convenience but has the feature of narrowing from b c to bd over a distance RS, which captures one significant geometric 
feature of the semicircular canal (Fig. lb). Neglecting hydrodynamic interactions with the wall, the sphere settles at 
the Stokes settling velocity, so that its axial position is given by s p (t) = U s t, where U s is given by Eq. (10). Viscous 
resistance to the settling sphere sets up a time-dependent pressure drop (Fig. IB3b ) on the cupula found by combining 
Eqs. (8), and CE), 

A , , 32a 3 Ann , . 

AP (t) = — (B2) 

3{(b c + b d )-(b c -b d )ta,nh[(U s t-s Q )/R5}} 2 ' 
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FIG. B3 (a) Trajectory of an otoconium falling through the straight, tapering canal described by Eq. IBlt . Heavy, filled dots 
denote the otoconium position every two seconds, (b) Cupular pressure AP c (t), displacement V c (t) and nystagmus E(t) as an 
otoconium falls along the centerline of the straight, tapering channel (whose profile is also indicated). The settling otoconium 
sets up a pressure field that is small when the otoconium is in the wide part of the channel (t < 3 seconds) and jumps by a 
factor b 2 c jb\ ~ 18 when the otoconium enters the narrow part of the channel. Correspondingly, the cupular volume displacement 
increases slowly while the force is small (w 3 seconds), then increases over a timescale r c to the large, constant force as the 
otoconium falls through the duct. The nystagmus response is similar, but increases over the longer time scale r v . Several 
seconds of 'latency' can be seen as the particle settles through the wide part of the channel, after which the force-amplifying 
effect of the narrowing canal sets in. Geometrical and physical parameters are given in Table I. 



As the particle settles through the narrowing portion of the channel, the pressure it exerts on the cupula increases 
by a factor b1/b 2 d ms 18, and asymptotically approaches a constant value AP C °° = 8a 3 Apg/3b d . The cupular volume 
displacement (Fig. IB3b l is then given by Eq. (7). After many time constants r c , the cupular volume displacement 
asymptotically approaches its limiting value 



AP c °°t c _ 8a 3 Apg 

— ~^KbJ 



V? = = « 57 P L > ( B3 ) 



where we have evaluated V c °° using the typical values reported in Table I. The nystagmus E(t) (Fig. IB3b ). calculated 
using Eq. (14), responds over the longer time scale r„, and asymptotically approaches its limiting value, 

\4 a poo 

^ = ,o h ±wkm ^ a28rad ' or 16 °- ( B4 ) 

This simulation is representative of basic hydrodynamic processes and corresponding cupular responses, but is never- 
theless unrealistic, as the real semicircular canal is curved and the forcing is truncated when the otoconium hits the 
wall. It does, however, clearly demonstrate the pressure-amplification effect of the tapering canal. 



APPENDIX C: Summary of results for a sphere settling in a circular cylinder 

This appendix summarizes several results for the hydrodynamic influence of the cylinder walls upon the motion of 
a small sedimenting particle. Our motivation for including this appendix is twofold: 1) to explain how hydrodynamic 
coupling between the wall and the particle can be treated in a reasonably straightforward yet accurate manner, and 
2) to demonstrate that this hydrodynamic coupling is largely unimportant in modelling the sedimentation of small 
particles in the semicircular canals. In regard to the second point, the walls only exert a significant influence upon the 
motion of a sedimenting particle when the particle is very close to the wall. In this limit, however, the fluid pressure on 
the cupula due to the sedimenting particle (given by Eq. (8)) is exceedingly small. Therefore, wall-particle interactions 
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FIG. C4 (a) Numerical factors for Eqs. ||UTJ giving corrections to the velocity of a sphere in an infinite, straight, circular 
cylinder which is forced radially (Wii) or axially (W33). Here f3 = r p /b is a measure of the sphere's fractional distance to the 
wall. Both Wn and W33 diverge like (1 — /3) _1 as the sphere nears the wall, which is the reason we have plotted (1 — (3)Wi\ 
and (1 — (3)W33. (From (Hirschfeld et al, 1984), Table 2.) (b) Settling velocity of a sphere parallel to the wall, in the limit 
where the sphere-wall gap is much smaller than the sphere radius (Eqs. IC2l and fC3t . 



are typically not significant for modelling cupular volume displacement, but are important when the trajectories of 
the particles are themselves of interest. 

A simple expression for the sedimentation velocity of a sphere in a geometry as complicated as that of the semi- 
circular canals is unavailable. However, because the radius of curvature R of the torus is much larger than the canal 
radius b, we approximate the canals as straight circular cylinders, which involves errors of order O (b/R). Even for 
this comparatively simple geometry, no uniformly valid expression is known for the motion of a sphere through a 
viscous fluid in a circular cylindrical container. Rather, numerical and asymptotic formulae have been derived in 
various limits. In this appendix, we summarize several results, which we use in our simulation of canalithiasis. 

We consider a small spherical particle of radius a, with externally applied force F g = F«z + F±r, located at radius 
r p in a circular cylinder of radius b. We introduce two dimensionless parameters, (3 = r p /b and k = a/6 to characterize 
the system. The parameter f3 indicates the dimensionless radial position of the particle, and k gives a dimensionless 
particle radius. 

In modelling the sedimentation of otoconia in the semicircular canals, we are concerned with the motion of particles 
that are small with respect to the size of the cylinder, and thus restrict our attention to the regime k <C 1. Within 
this limit, there are two primary regimes: 1) 1 -/? > k, where the distance from the particle to the wall is large 
compared to the particle size, and 2) 1 — (3 ~ k, where the particle is close to the wall. 

For the first regime, a method of reflections can be used to perturbatively calculate the influence of the wall upon 
the motion of the sphere. While results accurate to O (« 2 ) are available (Hirschfeld et al., 1984), we retain terms to 
O (k) only, with perpendicular and parallel velocities are given by 

C/±(r p »a) = -^(l- K W u ) and 17m (r p » a) = — L (1 - kW 33 ) , (CI) 
07r/ia 1 bnfia 

where W11 and W33 are dimensionless functions (Fig. IC4J) . 

A different behavior is found when the sphere is close to the wall, and the dimensionless gap e — (b — r„)/a is 
small. The leading-order approximation to the settling speed of the sphere can be found by treating the cylinder wall 
as locally planar, in which limit asymptotic formulae have been determined using a lubrication approach (Cooley & 
O'Neill, 1968; Goldman et al. , 1967; O'Neill & Stewartson, 1967). The formula given by Goldman et al. contains an 
error, and should read 



" 67T/iO 



1.91-21n(e) 



1.59-3.191n(e) + {ln(e)} 2 



(C2) 



which very nearly agrees with the results obtained in (Cooley & O'Neill, 1968; O'Neill & Stewartson, 1967), which 
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give 



67T/iO 



1.85-21n(e) 



1.52 - 3.15 ln(e) + {m(e)} 2 



(C3) 



The velocity parallel to the wall thus decays logarithmically slowly. 

Motion perpendicular to a planar wall due to a force F± can be found using the lubrication approximation (e.g. 
Batchelor, 1967, p. 228, problem 1) in the limit e <C 1, giving 



F ± e 



(C4) 



so that a sphere exponentially approaches the wall with time constant t± = 6nfia 2 /F±. 
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